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We analyze statistical properties of complex eigenvalues of 
random matrices A close to unitary. Such matrices appear 
naturally when considering quantized chaotic maps within a 
general theory of open linear stationary systems with discrete 
time. Deviation from unitarity are characterized by rank M 
and eigenvalues Ti,i = 1, M of the matrix T = 1 — A' A. 
For the case M = 1 we solve the problem completely by de- 
riving the joint probability density of eigenvalues and calcu- 
lating all n— point correlation functions. For a general case 
we present the correlation function of secular determinants. 

PACS numbers: 03.65 Nk, 05.45 Mt 



The theory of wave scattering can be looked at as an in- 
tegral part of the general theory of linear dynamic open 
systems in terms of the input-output approach. These 
ideas and relations were developed in system theory and 
engineering mathematics many years ago, see papers 
and references therein. Unfortunately, that develop- 
ment went almost unnoticed by the majority of physicists 
working in the theory of chaotic quantum scattering and 
related phenomena, see and references therein. For 
this reason I feel it could be useful to recall some basic 
facts of the input-output approach in such a context. 

An Open Linear System is characterized by three 
Hilbert spaces: the space Eq of internal states \& € Eq 
and two spaces E± of incoming (-) and outgoing (+) sig- 
nals or waves also called input and output spaces, made 
of vectors <j>± E E±. Acting in these three spaces are 
four operators, or matrices: a) the so-called fundamental 
operator A which maps any vector from internal space 
Eq onto some vector from the same space Eq b) two op- 
erators Wi 2) with W\ mapping incoming states onto an 
internal state and W<i mapping internal states onto out- 
going states and c) an operator Sq acting from E_ to 
E+. 

We will be interested in describing the dynamics 
of an internal state with time t provided we know the 
state at initial instant t = and the system is sub- 
ject to a given input signal 0_ (t) . In what follows we 
consider only the case of the so-called stationary (or 
time-invariant) systems when the operators are assumed 
to be time-independent. Let us begin with the case of 
continous-timc description. The requirements of linear- 
ity, causality and stationarity lead to a system of two 
dynamical equations: 



AAt(t) + Wif-(t) 



(1) 



Interpretation of these equations depends on the nature 
of state vector '3/ as well as of vectors <fi± and is differ- 
ent in different applications. In the context of quantum 
mechanics one relates the scalar product with the 
probability to find a particle inside the " inner" region at 
time t, whereas <$>± ( t > ± stays for probability currents flow- 
ing in and out of the region of internal states (the number 
of particles coming or leaving the inner domain per unit 
time) . The condition of particle conservation then reads 
as: 



dt 



(2) 



It is easy to verify that Eq. (j2|) is compatible with the 
dynamics Eq.(|l|) only provided the operators satisfy the 
following relations: 



- A = iWW^ 



SISq = i andTF f = W 2 = -S W} 



which shows, in particular, that A can be written as A = 
H - \WW\ with a Hermitian H = Hi 

The meaning of H is transparent: it governs the evo- 
lution i-§{^ — H^{t) of an inner state \& when the 
coupling W between the inner space and input/output 
spaces is absent. As such, it is just the Hamiltonian 
describing the closed inner region. The fundamental op- 
erator A then has a natural interpretation of the effective 
non-self adjoint Hamiltionian describing the decay of the 
probability from the inner region at zero input signal: 
(j>-(t) — for any t > 0. If, however, the input signal is 
given in the Fourier-domain by <fi-(u>), the output signal 
is related to it by: 



4>+{ui) = S{u>)S 4>-{lo) 



iW* 1 - 



1 



-W 



ujl-A 
(3) 



where we assumed 9(t = 0) = 0. The unitary matrix 
S(w) is known in the mathematical literature as the char- 
acteristic matrix-function of the non-Hermitian operator 
A. In the present context it is just the scattering matrix 
whose unitarity is guaranteed by the conservation law 
Eq®. 

The contact with the theory of chaotic scattering is 
now apparent: the expression Eq.(^|) was frequently used 
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in the physical literature as a starting point for ex- 
tracting universal properties of the scattering matrix for 
a quantum chaotic systems within the so-called random 
matrix approach. The main idea underlying such an ap- 
proach is to replace the actual Hamiltonian H by a large 
random matrix and to calculate the ensuing statistics of 
the scattering matrix. The physical arguments in favor of 
such a replacement can be found in the cited literature. 

In particular, most recently the statistical properties 
of complex eigenvalues of the operator A as well as re- 
lated quantities were studied in much detail |^-|^] . Those 
eigenvalues are poles of the scattering matrix and have 
the physical interpretation as resonances - long-lived in- 
termediate states to which discrete energy levels of the 
closed systems are transformed due to coupling to con- 
tinua. 

In the theory presented above the time t was a con- 
tinous parameter. On the other hand, a very useful in- 
strument in the analysis of classical Hamiltonian systems 
with chaotic dynamics are the so-called area-preserving 
chaotic maps [|9| . They appear naturally cither as a map- 
ping of the Poincare section onto itself, or as a result of 
stroboscopic description of Hamiltonians which are peri- 
odic function of time. Their quantum mechanical ana- 
logues are unitary operators which act on Hilbert spaces 
of finite large dimension N . They are often referred to as 
evolution, scattering or Floquet operators, depending on 
the physical context where they are used. Their eigenval- 
ues consist of N points on the unit circle (eigenphases). 
Numerical studies of various classically chaotic systems 
suggest that the eigenphases conform statistically quite 
accurately the results obtained for unitary random ma- 
trices (Dyson circular ensembles). 

Let us now imagine that a system represented by a 
chaotic map ("inner world") is embedded in a larger 
physical system (" outer world" ) in such a way that it 
describes particles which can come inside the region of 
chaotic motion and leave it after some time. Models of 
such type appeared, for example, in [f[o| where a kind of 
scattering theory for "open quantum maps" was devel- 
oped based on a variant of Lipmann-Schwinger equation. 

On the other hand, in the general system theory [|]]|] 
dynamical systems with discrete time are considered as 
frequently as those with continous time. For linear sys- 
tems a "stroboscopic" dynamics is just a linear map 
(<j)-{n); 4'(n)) — > (0+(n); ^(n + 1)) which can be gen- 
erally written as: 

¥(n + l)\ =ir f *(n) \ # (A Wi 
d>+(n) J \4>-{n) J ' \W 2 S 

(4) 

Again, we would like to consider a conservative system, 
and the discrete-time analogue of Eq.(|^) is: 

\E ,t (n+l)*(n+l)-* t (n)*(n) = <f^_{n)(f>-(n)—(f)\. (n)<j>+(n) 

which amounts to unitarity of the matrix V in Eq.(|]). 



In view of such a unitarity V of the type entering Eq. 
can always be parametrized as (cf. [El)): 



V 



ui 
vi 



VT 



-v^i-ftf) 



U2 

v 2 



(•5) 



where the matrices u± t 2 and v±.2 are unitary and f is 
a rectangular N x M diagonal matrix with the entries 
m = SijTi, l<i<N,l<j<M <; tj < 1. 

Actually, it is frequently convenient to redefine 
input, output and internal state as: 4>-{n) — > 
^2^- 4>+ ~ * vi4'+( n ) an d Vl'(n) — > U2^(n) which 
amounts just to choosing an appropriate basis in the cor- 
responding spaces. The transformations bring the matrix 
V to somewhat simplier form: 



V 



iVT 



rt 



-VT 



(6) 



where u = u\u\. Such a form suggests clear interpreta- 
tion of the constituents of the model. Indeed, for f = 
the dynamics of the system amounts to: ^{n + 1) = 
iity(n). We therefore identify u as a unitary evolution op- 
erator of the " closed" inner state domain decoupled both 
from input and output spaces. Correspondingly, f ^ 
just provides a coupling that makes the system open and 
converts the fundamental operator A — uV 1 — rrt to a 
contraction: 1 — A' A = tt^ > 0. As a result, the equa- 
tion \&(n + 1) = A^(n) describes an irreversible decay of 
any initial state ^(0) for zero input <fi-(n) — 0, whereas 
for a nonzero input and ^(0) = the Fourier-transforms 
(f)±(uj) — J2^=o e muJ (t>±(n) are related by a unitary scat- 
tering matrix S(u>) given by: 



S(lu) = -VI 



T'T + T 



1 



- A 



(7) 



Assuming further that the motion outside the inner 
region is regular, we should be able to describe generic 
features of open quantized chaotic maps choosing the 
matrix u to be a member of a Dyson circular ensem- 
ble. Then, averaging S(w) in Eq.(0) over u one finds: 

Comparing this result with PI 
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we 

see that M eigenvalues < T$ < 1 of the M x M matrix 
T = f'T play a role of the so-called transmission coeffi- 
cients and describe a particular way the chaotic region is 
coupled to the outer world. 

In fact, this line of reasoning is motivated by recent pa- 
pers IT^Il^ The authors of [jTJj considered the Floquet 
description of a Bloch particle in a constant force and pe- 
riodic driving. After some approximations the evolution 
of the system is described by a mapping: c, i+ i = Pc n , 
where the unitary Floquet operator F = SU is the prod- 
uct of a unitary "M-shift" S : Ski = &i,k-M, l,k = 
— oo,oo and a unitary matrix U. The latter is effec- 
tively of the form U = diag(di, u, d^), where di t 2 are 
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(semi) infinite diagonal matrices and u can be taken from 
the ensemble of random N x N unitary matrices. 

One can check that such a dynamics can be easily 
brought to the standard Eqs.(|],^|) with the fundamen- 
tal operator being N x N random matrix of the form 
A = vl — Wru, and all M diagonal elements of N x M 
matrix r are equal to unity. Actually, the original paper 
fl2"f employed a slightly different but equivalent construc- 
tion dealing with an "enlarged" internal space of the di- 
mension N + M. We prefer to follow the general scheme 
because of its conceptual clarity. 

Direct inspection immediately shows that the non- 
vanishing eigenvalues of the fundamental operator A as 
above coincide with those of (N — M) x (N — M) sub- 
block of the random unitary matrix u. Complex eigen- 
values of such "truncations" of random unitary matrices 
were studied in much detail by the authors of a recent 
insightful paper |l3[ ] . They managed to study eigenvalue 
correlations analytically for arbitrary N, M . In partic- 
ular, they found that in the limit N — * oo for fixed 
M these correlation functions practically coincide JlJ] 
with those obtained earlier for operators of the form 
A = H — ^WW^ occuring in the theory of open systems 
with continous-time dynamics. 

Such a remarkable universality, though not completely 
unexpected, deserves to be studied in more detail. In 
fact, truncated unitary matrices represent only a par- 
ticular case of random contractions A. Actually, some 
statistical properties of general subunitary matrices were 
under investigation recently as a model of scattering ma- 
trix for systems with absorption, see |l4j] . 

The main goal of the present paper is to add to our 
knowledge on specta of random contractions for a given 
deviation from unitarity. 

The ensemble of general N x N random contractions 
A = uy/1 — tt t can be described by the following prob- 
ability measure in the matrix space: 



T(A)dA tx S(A^A - G)dA , G = 1 - ff 1 



(8) 



rt = 



where dA = J]^ dA^dA*^. The N x N matrix fft 

1 — G > is natural to call the deviation matrix. It has 
M nonzero eigenvalues coinciding with the transmission 
coefficients T a introduced above. The particular choice 
Tj <M — 1, T;>Af = corresponds to the case considered 
in |l3f|. In what follows we assume all Tj < 1. 

Our first step is, following [||Jl^], introduce the Schur 
decomposition A — U(Z + R)ffi of the matrix A in terms 
of a unitary U, diagonal matrix of the eigenvalues Z and 
a lower triangular R. One can satisfy oneself, that the 
eigenvalues zi, zjv are generically not degenerate, pro- 
vided all Tj < 1. Then, the measure written in terms 
of new variables is given by dA = \A({z})\ 2 dRdZd^i(U), 
where the first factor is just the Vandermonde determi- 
nant of eigenvalues Zj and d/i(U) is the invariant measure 
on the unitary group. The joint probability density of 
complex eigenvalues is then given by: 



P({z})(x|A({z})| 2 (9) 
x / dn(U)dR 5 {{Z + R)(Z + - U^G&) 

The integration over R can be performed with some ma- 
nipulations using its triangularity (some useful hints can 
be found in fll3l). As the result, we arrive at: 



V({z}) cx |A({z})| ; 



(10) 



det 



WGU 



(»,i)=i,...,j, 



The remaining integration over the unitary group poses 
a serious problem. We found no way to overcome the 
difficulties in a general case. However, when the rank 
M of the deviation matrix tt* is unity, we managed to 
perform the integral by methods of Q and arrived at a 
very simple expression: 



^({zDcxT^IAaz})! 2 ^ 



1 — T 



(ii) 



provided < |z/| < 1 for all eigenvalues, and zero other- 
wise. Here < T < 1 is the only non-zero eigenvalue of 
the deviation matrix. The Eq. ( pT| ) can be used to extract 
all n-point correlation functions: 

i?„(zi, z n ) = /A Nl f d 2 z n+1 ...d 2 z N V({z}) (12) 



(N-n)\ 

To achieve this it is convenient to use the Mellin trans- 
form with respect to the variable C = 1 — T: 



Rn(s; {z} n ) 



dCC-MCi-O^^Wn)] (13) 



It easy to notice that such a transform brings V({z}) 
to the form suitable for exploitation of the orthogonal 
polynomial method [ fl5[ . The corresponding polynomials 
are Pk{z) = y/Jk + s)z n orthonormal with respect to the 
weight f(z) — |z| 2 ( s_1 ^ inside the unit circle \z\ < 1. 
Following the standard route we find: 



Rn(s; {z} n ) OC 



det [K(zi, 



z j)l\(i,j)=l, 



s(s + l)...(s + N -1) 



where the kernel is 

K(z 1 ,z 2 ) = (f(z 1 )f(z 2 )) 



N-l 

^ X>*(*r)p*te) 

fc=0 



(14) 



(15) 



\x\ s 1 (s4>(x) + x^-(t){x)^ \ X = Z{Z2 



and 4>(x) — (x — l)/( x — !)• Thus, the expression 
Eq.(O) can be rewritten as: 



R n (s; {z} n ) cx 
qo({z}n) 

Qn({z}n) 



nr=ii^i 



s(s + l)...(s + N-l)^ 



l Qi({z} n ); 



det 
det 



X -^, ( t>i. x )\x=(z'z j ) 
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and can be easily Mellin-inverted yielding finally the orig- 
inal correlation functions in the following form: 

a 

Rn{{*}n)\M<i « T 1 - N 9(T -l + a)J2 ffl({*}n) (16) 



x I —a 

, da 



1L_1_T 
a V a 



1=0 

N-l 



a =IYU i 2 



where 9(x) = 1 for x > and zero otherwise. This 
equation is exact for arbitrary TV. Let us now investigate 
the limit N ^ n and use: 



da 



1 r 



a V a 



JV-l 



a — £/ a \ a 



AT-1 



which allows one to rewrite the correlation function as: 

JV-l 

0(T-l + a) (17) 



^({^txT 1 ^- 1 
a 



1 — T 



/JV(l-T) . . rf.A. 

x det 7—^^ x + ^r^ 1 

j,j=i,...,n \a — 1 + 1 dx J 

Further simplifications occur after taking into account 
that eigenvalues z% are expected to concentrate typically 
at distances of order of 1/N from the unit circle. Then 
it is natural to introduce new variables yi,fa according 
to Zj = (1 — yi/N)e l ^ i and consider to be of the order 
of unity when N — > oo. As to the phases fa, we expect 
their typical separation scaling as: fa — 4>j = 0(1/N). 
Now it is straightforward to perform the limit N — > oo 
explicitly and bring Eq.dTTl) to the final form: 



R n ({z}„) oc e 



det 



d\(X + g)e 



-■k\s, 



with g = 2/T — 1 and 5y = N(fa — <j)j) — i(yi + yj). 
The expression above coincides in every detail with that 
obtained in for random matrices with rank-one de- 
viation from Hermiticity provided one remembers that 
mean linear density of phases fa along the unit circle is 
v = l/(27r). This completes the proof of universality for 
rank-one deviations. 

Being so far unable to evaluate the spectral correlation 
function for an aritrary A, we nevertheless succeeded in 
calculating closely related but simpler object, namely the 
correlation function of secular determinants: 



I(z u z 2 ) = ^(dct (zil - A\ det (41 - A^ 



(18) 



where the angular brackets stand for the averaging over 
the probability density in Eq.(||). Such a correlation 
function for unitary A corresponding to quantum chaotic 
maps attracted much attention recently ]9|,|l6[| . For a non- 
Hermitian A similar objects were studied in [f|. 

Relegating details to a more extended publication, we 
just present the final result in terms of eigenvalues gi = 
1 — Ti of the matrix G: 



I(zi,z 2 ) 



N 



E( z ^*) fe 



fe=0 



N 



E 



(19) 



9i\9ii ■■•9ik 



Kil<i 2 <—<ik<N 



For random unitary matrices all Tj = and the expres- 
sion above reduces to I{zi,Zi) = X)^=o( z i z 2)' c m agree- 
ment with p6[ . 
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